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ABSTRACT 

We extend the methods of Yennie, Frautschi and Suura (YFS) to compute, via Monte Carlo 
methods, the effects of multiple gluon emission in the processes q + { (f' — > q" + ( q ) "' + n(G), 
where G is a soft gluon. We show explicitly that the infrared singularities in the respective 
simulations are canceled to all orders in a s . Some discussion of this result from the standpoint 
of confinement is given. More importantly, we present, for the first time ever, sample numerical 
Monte Carlo data on multiple soft gluon emission in the rigorously extended YFS framework. 
We find that such soft gluon effects must be taken into account for precise SSC/LHC physics 
simulations. 
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1 Introduction 



While the SSC has been cancelled and is being re-evaluated, the LHC planning stages continue 
to gain momentum, so that it still becomes more and more necessary to prepare for the physics 
exploration which the LHC, at the least, will provide. The primary issue in this exploration 
is the comparison between predictions from the standard SU2L ® TJ\ (g> SU^ theory and its 
possible extensions, and what will be observed by the GEM and SDC collaborations, with the 
discovery of the Higgs particle or whatever it represents as of course a primary goal of such 
comparisons. Thus, it is important to know the theoretical predictions as accurately as is 
necessary to exploit the expected detector performances of GEM and SDC to the fullest extent. 
In particular, higher order radiative corrections to these predictions (due to multiple photon 
and multiple gluon effects) are in fact essential to obtain the proper precision on the signal and 
background processes in the respective SSC/LHC environment. In Refs. [|IJ, we have developed 
Yennie-Frautschi-Suura |2| (YFS) multiphoton Monte Carlo event generators for calculating 
nj effects in SSC/LHC processes. In what follows, we now develop the extension of our own 
YFS Monte Carlo-based higher order radiative correction methods to multiple gluon effects in 
SSC/LHC processes. Such a calculation of n(G) effects has not appeared elsewhere, where we 
use G to denote a gluon. 

More precisely, we want to use the fact that, by the uncertainty principle, infinite-wavelength 
gluons should not affect the motion of quarks and gluons inside a proton, whose radius is ~ 1 
fm. Thus, we have recently realized || this physical requirement in perturbation theory by 
showing that in our prototypical processes q + — ► q" + { q ) "' + (G), the infrared singularities 
cancel at 0(a s ), just as they cancel in the analogous QED process f + f — > f + f' + J- We then 
are able to define the QCD analogues of the famous YFS infrared functions || B and B which 
describe the virtual and real infrared singularities in QED processes to all orders in a. These 
functions, -Bqcd and -Bqcd, then describe the infrared singularities in QCD in the processes 
q+ ( q ) — > q" + ( q ) "' + n{G) to all orders in a$- The sum £?qcd + -Bqcd is infrared finite and allows 
us to extend our YFS Monte Carlo methods, widely used in SLC/LEP physics for higher order 
QED radiative effects, to the analogous Monte Carlo methods in QCD processes for SSC/LHC 
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physics objectives. We call the resulting programs SSCYFSG and SSCBHLG, for they are the 
QCD extensions of the QED w-f Monte Carlos SSCYFS2 and SSCBHL. Here we recall that 
SSCYFS2 was already described in Ref. |]J and SSCBHL is the analogous SSC extension of the 
YFS Monte Carlo BHLUMI 2.01 g for the SLC/LEP luminosity process e+e" -> e+e~ + nj. 

In what follows, we shall present sample Monte Carlo data for the process q + — > 
q" + { q ) "' + n(G) from the event generator SSCBHLG. Analogous results from SSCYFSG will 
appear elsewhere [H. The latter Monte Carlo only simulates initial-state gluon radiation, 
whereas SSCBHLG simulates initial-state and final-state n-gluon radiation, as well as the 
respective initial-final-state interference effects. Further, we emphasize that (unlike BHLUMI 
2.01) SSCBHLG is not restricted to small scattering angles. Note finally that SSCYFSG is an 
initial-state restriction of SSCBHLG and hence is useful for cross-checking our work. 

Our work is organized as follows. In the next Section, we derive the extension of the YFS 
exponentiation to n-gluon emission and the attendant extension of SSCYFS2 and SSCBHL. In 
Section III, we present sample Monte Carlo data for n-gluon emission for the SDC and GEM 
acceptances. Section IV contains some concluding remarks. 

2 Infrared Singularity Cancellation in QCD 

In this Section, we derive the analogue of the YFS infrared functions B and B for QCD, 
Bqcb and -Bqcd, for the process q + — ► q" + ( g } + (G). We will focus on the case of 
most interest to us, in which a gluon is exchanged in the i-channel and we require that the 
exchanged gluon carries large momentum transfer so that the outgoing quark (anti-quark) is 
in the 1 77 1 < 2.8 acceptance of the SDC (and GEM). This situation is depicted in Fig. |l|. We 
remark that if the relevant exchange is a color singlet exchange (such as W or Z° exchange), 
we can obtain the corresponding formulas for -Bqcd and -Bqcd by dropping the terms in the 
color-exchange case which arise from non-commutativity of the respective QCD A a matrices in 
the relevant quark representation. 

Specifically, following the kinematics in Fig. |I]a, we note that, for the SDC and GEM 
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acceptance, we always have — Q 2 = — (^i — (fc) 2 ^ Aqcd> so that perturbative QCD methods are 
applicable. What we wish to do is to extract the infrared-singular part of the 0(a s ) amplitude 
in Fig. [l|b, in complete analogy with the YFS extraction of the infrared-singular portion 
of the analogous amplitude in QED via the famous YFS virtual infrared function B. To this 
end, we first note that, by explicit calculation, the graphs (v)-(vii) are not infrared-divergent 
(we compute the gluon vacuum polarization to order a s as well), so that the infrared-singular 
part of Fig. |l]b is given by just the same graphs as the infrared-singular part of the analogous 
process in QED. We follow the YFS methods in Ref. p| and find the infrared-singular part of 
Fig. |l]b to be given by 



•M Vj ir = a s B QC T)MsB + a s B QCD MoB : 

where M.%b is the respective Born amplitude in Fig. |l]a and where 
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with C F = 4/3 =quadratic Casimir invariant of the quark color representation, A equal to the 
standard infrared regulator mass (Ref. 0) and 
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Thus, we see that the non-Abelian nature of QCD causes ACj ^ 0. 

The additional term in (|]), which arises from color algebra, consists of what we call the 
"singlet" exchange Born amplitude which is obtained from Fig. |l|a by substituting A for the 
color matrices A a , where A = I/v6. The corresponding value of -Bqcd in ([[]) is obtained from 
(0) by setting Cp = and AC t = in (0) and by setting 

AC U = AC S = { qq \ inC ° ming (4) 

I Cf, qq incoming 

in (0). We note that, to order a s , the "singlet" amplitude does not contribute to our cross 
section for the process in Fig. |]. 

Turning now to the soft real emission process in Fig. [l]b, we again follow the methods of 
YFS in Ref. and extract the infrared-divergent part of the real gluon bremsstrahlung as 

dafaft = da (2a s B Q cD) (5) 

where the real QCD infrared function is 

2a s B QC D = / -j— SqcdO) (6) 

with 
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where K max corresponds to the relevant gluon jet detector resolution energy (~ 3 GeV at 



the SSC/LHC) and fco = yk 2 + \ 2 so that we regulate the infrared singularities in (0) in 
complete analogy with our gluon mass regulator in (0). (Note that graph (v) in Fig. |l|c is 
not infrared-divergent.) Here, dao is the respective Born differential cross section. The results 
(|H) to (0) represent the complete infrared (IR) structure of QCD for the cross section for 



q + — ¥ q" + ( g )W + (G) to order a s . The result (]?]) is consistent with the pioneering work 
of Berends and Giele in Ref. J7|. Following the YFS theory, we see that the criterion for the 
cancellation of the 0(a s ) and IR singularities in the cross section for our scattering process in 
Fig. is that 

SUM IR (QCD) = 2a s B QC B + 2a s ReB QC v (8) 

should be independent of A. Using the well-known results [@],|| for the integrals in (fj) and 
(0) (notice that these are just the same integrals that arise in the YFS analysis of QED with 
different color-based weights) we see that, indeed, A does cancel out of We are left with 
the fundamental result (for, e.g., m q = m q i = m) 



SUM IR (QCD) = -JL £ (-iyW(C F B tot (A) + AC A B' tot (A)) (9) 

71 A={s,t,u,s',t',u'} 



where 

B tot {A) = \og(2K max /J\A\) 2 (\n(\A\/m 2 ) - 1) + \ ln(|A|/m 2 ) - 1 - tt 2 /6 + 9(A)n 2 /2 



and 



Bf M (A) = log(2A moJ ,/y1!j) 2 + iln(|^|/m 2 ) - i 2 /6 + 9(A)^/2 (10) 



The results (0) and (|T0|) are then the fundamental results of our analysis. 

Indeed, repeating the arguments in Ref. @ and/or using the factorization results in Ref. 0, 
we see that by virtue of (^|), we may compute the hard gluon radiation residuals at 0(a s ) as 

P Q = da x - loov - 2a s ReB QCD d* (12) 

and 
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in the exponentiated formula 
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Here, we use the virtual 1-loop correction 
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. It allows (3q to be identified as da®}? s; 



2a s {Bqcd + B QC D)da , where da® o y a> 



da l-loop + 



da f t , for example. The above expressions are for incoming gg and need to be modified ac- 
cordingly for qq interactions. The hard gluon residual (3\ makes a vanishing contribution in the 
soft gluon limit so that it will be presented elsewhere ||; for, here we focus on the soft gluon 
effects in ( |I4D so that we will not need j3\ in the current discussion. 



The results 



and (061 



If 



are now realized in fll4|) via Monte Carlo methods by 
replacing the QED forms of j3 n , SUMir as described in Refs. [0 for our SSCYFS2 and SSCBHL 
event generators by the respective QCD forms. The resulting event generators are SSCYFSG 
and SSCBHLG, where in the former, only only gluon radiation from the initial state is treated. 
Sample Monte Carlo data from SSCBHLG are illustrated in the next section. Similar data for 
SSCYFSG will appear elsewhere Ifl. 
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We should stress that two limits of the QCD coupling constant are needed in (|T2|), ( |T3"D and 
(PI). One is the perturbative — q 2 ^> Aq CD regime in the hard gluon residuals j3 n in (|T3|). This 
regime is well-known and we use the standard type of formula 



a s (a) = — == — (19) 

(33 - 2n f ) ln(/x 2 /(A^ s ) 2 ) 



with rif =number of open flavors and A^f s the respective value of Aqcd : we take A% IS ~ 
.238 GeV for definiteness. 

The second regime occurs in SUMjr in @, where the k 2 — > limit is relevant for "on-shell" 
gluons of 4-momentum k. (We follow Tarrah [jnj and use the concept of on-shell quarks and 
gluons in perturbation theory at large momentum transfer.) Recently ||11|| , it has been pointed 



out by Mattingly and Stevenson that this limit of a s /ir exists as an infrared fixed point and its 
limiting value is ~ .263. We use this value of (a s /ii) in SUM IR in @ and ([U|). Because of the 
fixed point behavior, the value one uses for the k 2 — > coupling of gluons is not sensitive to the 
precise value A4 , for example. Further, since we only work to 0(a s ) in the hard interactions, 
there is no contradiction between our use of A MS values in ( |T9"D and the renormalization scheme 



of Mattingly and Stevenson in Ref. |TT 



We should also note that our light quark masses are always those determined by Leutwyler 
and Gaisser in Ref. [0]: m u (l GeV) ~ 5.1 x 1(T 3 GeV, m d (l GeV) ~ 8.9 x 1(T 3 GeV, 



m s (l GeV) ~ .175 GeV, m c (l GeV) ~ 1.3 GeV, mb(m b ) ~ 4.5 GeV; for m t , we take m t (m t ) ~ 

we 



164 GeV. The running of these masses is readily incorporated into our calculations |T2| 
ignore this running in the current paper, since it does not affect our results at the level of 
accuracy of interest to us here. 

With these explanatory remarks in mind, we now turn to explicit Monte Carlo data illus- 
trations for multiple gluon effects in q + — > q" + ( g )W + n(G) at SSC/LHC energies in the 
ATLAS-CMS (SDC-GEM) acceptances. 
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3 Results 



In this Section, we illustrate our multiple-gluon Monte Carlo event generators in the sample 
case ofw + w— > u + u + n(G), where we require the out-going u quarks to satisfy \r]\ < 2.8 
and E > 58 GeV, for definiteness, at y/s = 15.4 TeV. (We comment about the SSC case yfs = 
40 TeV as well.) For completeness, we will show explicit Monte Carlo data for SSCBHLG, since 
in that generator initial, initial-final interference, and final state n(G) radiation is simulated. 

Our results are shown in Figs. |2], |3| and where we show respectively the gluon multiplicity, 
the distribution of v = (s — s') / s, where s' is the squared invariant mass of the final quark-quark 
system, and the total gluon transverse momentum in TeV units. What we see is that there is a 
pronounced effect from the multi-gluon radiation at = 15.4 TeV: the average value of the 
number of radiated gluons with energy E G > 3 GeV is, e.g., at ^/s = 15.4 TeV, 

(n G ) = 28.5 ±5.5. (20) 
The average value of v is, for yfs = 15.4 TeV, 

(v) = 0.92 ± 0.14. (21) 
The average value of the total transverse momentum in gluons is 

n 

(j>±,tot) = (lE^D = (°- 19± °- 23 ) TeV - ( 22 ) 

i=l 

This means that a substantial amount of incoming and outgoing energy is radiated into gluons 
(due to the usual collinearity of radiation with its source, most of the final state radiation 
would be a part of the typical jet associated with its parent quark). Entirely similar results 
hold for the SSC y/s = 40 TeV case. Evidently, any realistic treatment of QCD corrections 
to corrections to LHC/SSC processes must analyze the full event -by-event n-gluon effects as 
they interact with the detector efficiencies and cuts. Our SSCBHLG event generator provides 
the first amplitude-based exponentiated Monte Carlo realization of such n-gluon effects and 
we hope to participate in a study of their interplay with detector effects elsewhere . 

An important consequence of our work is a prediction for the effect on the overall normal- 
ization of an LHC/SSC physics process due to multiple-gluon radiation, in an unambiguous 
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way, since all IR-singular effects have canceled out of our calculation. For example, if we sim- 
ply compute the 0(a s ) soft cross section, E G < 3 GeV, we get (for u + u — > u + u + (G), 
= 40 TeV) that 

&/&Born = 1 + &vir + 5 rea l = —22.1, (23) 

whereas if we compute the cross section for u+u — > u+u+n{G) and require v < 3 GeV/7.7 TeV= 
.00039, we get 

o /o Born = 1.2 x 10" 11 . (24) 

Hence, we see that multi-gluon radiation is crucial to getting the proper normalization of the 
cross section. 

We emphasize that the results for = 40 TeV give similar conclusions || to those in 
Figs. |2|-44|. Entirely analogous results hold for the pure initial-state multigluon event generator 
SSCYFSG, except only two lines radiate, so that (n G ) is reduced by a factor ~ 2, with the 
appropriate corresponding reductions in (v) and (pj_ t tot)- Such results are effectively an inte- 
gration over the events used in Figs. |2|-44|, wherein one integrates over all gluons within some 
cones close to the outgoing fermions. We do not show a separate set of plots for SSCYFSG here 
but we will present a more detailed discussion of the respective initial-state n-gluon radiation 
elsewhere ||. Here, we are focusing on the complete gluon radiation problem. 



4 Conclusions 



In this paper, we have shown that the IR singularities in fermion-(anti-)fermion scattering 
at LHC/SSC energies in QCD cancel at order a s , permitting an immediate extension of YFS 
QED exponentiation methods to such processes. The YFS Monte Carlo methods invented by 
two of us |T4| (S.J. and B.F.L.W.) for QED radiation in Z° physics can then be extended to 
multi-gluon Monte Carlo event generators with well-defined IR behavior to all orders in a s . 
Such a Monte Carlo realization of amplitude-based n-gluon effects opens the way to a new era 
in higher-order QCD corrections to hard processes such as those of interest to the LHC/SSC 
physics program. 
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More specifically, we have realized our exponentiated, n-gluon effects via the Monte Carlo 
event generators SSCBHLG and SSCYFSG. Since the latter amounts to an initial-state re- 
striction of the former, we have presented results from the former generator in this paper. We 
find, for example, that the average number of radiated gluons in u + u — > u + u + n(G) at 
y/s = 15 A TeV is 28.5 ± 5.5 and that (v) = 0.92 ± 0.14. Further, the cross-section normal- 
ization is strongly affected by our n-gluon effects. Hence, only by taking these effects into 
account, as they interact with detector simulation, for example, can a truly detailed picture of 
the SSC/LHC physics platform be obtained. We hope to participate in such an investigation 
elsewhere. 

In summary, we have presented, for the first time, the realistic event-by-event Monte Carlo 
realization of multi-gluon effects in hard QCD processes at SSC/LHC energies in which IR 
singularities are canceled to all orders in a s . A new approach to LHC/SSC physics has thus 
been created. We look forward with excitement to its many applications! 
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(i) (11) (ill) 




(c) 

Figure 1: The process q + { q )l — > q + ( q }/ + (G) to 0(a s ): (a) Born approximation; (b) 0(a s ) 
virtual correction; (c) 0(a s ) bremsstrahlung process. Here, f.s.e. represents the fermion self- 
energies. 
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Figure 2: Gluon multiplicity distribution mu + u^u + u + n(G) at y/s = 15.4 TeV. 
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Figure 3: f-distribution in u + u^u + u + n(G) at ^fs = 15.4 TeV. 
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Figure 4: Total transverse momentum of gluons in u + u—*u + u + n(G) at y/s = 15.4 TeV. 
The units of the horizontal axis are TeV. 
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